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The  correct  design  of  underwater  camera  systems  for  viewing  underwater  objects  is 
vitally  important  if  the  performance  of  these  underwater  imaging  systems  is  to  be  maximized. 
Towards  this  goal,  the  Woods  Hole  Oceanographic  Institution  has  developed  a  system  of 
computer  programs  which  allows  the  underwater  lighting  system  designer  to  explore  the 
imaging  system  performance  that  results  from  the  manipulation  of  beam  patterns,  geometry 
of  cameras  and  light  sources,  and  changes  in  the  environment. 

The  computer  simulation  of  underwater  image  system  performance  has  been  found  to  be 
a  valuable  tool  for  several  reasons.  Most  importantly,  the  performance  of  underwater  lighting 
systems  cannot  be  easily  predicted  from  terrestial  experience  because  of  the  intense  scattering 
of  the  oceanic  medium.  Secondly,  the  cost  of  implementing  and  running  computer  programs 
to  simulate  underwater  camera  light  viewing  is  at  greatly  reduced  expense  to  experimentation 
in  the  real  world.  Finally,  the  inherent  flexibility  in  using  a  computer  for  modeling  allows  the 
user  to  build  up  a  base  of  experience  which  can  then  be  used  for  heuristic  system  design. 

In  this  report,  we  will  consider  the  results  of  a  systematic  study  that  was  performed 
in  order  to  quantify  and  refine  the  performance  of  an  underwater  imaging  system.  The 
camera  and  lighting  system  of  the  mine  neutralization  system  (MNS)  was  subject  to  an 
intensive  computer  study  of  over  500  simulations  in  order  to  characterize  the  existing  system 
performance  and  to  determine  the  scope  of  both  simple  and  more  complex  changes  that  could 
be  made  in  order  to  optimize  the  performance  of  the  imaging  system. 

Figure  1  is  a  page  from  a  document  which  details  the  MNS  performance.  As  documented,  it 
is  evident  that  the  mission  requirement  of  the  vehicle  is  to  neutralize  either  bottom  or  moored 
mines,  once  they  have  been  located  by  some  other  device.  In  the  bottom  neutralization  mode, 
the  MNS  deploys  a  small  bomblet  within  a  given  radius  of  the  mine  which  is  then  detonated 


$ 
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after  the  vehicle  has  swum  a  safe  distance  away.  The  primary  requirement  here  is  that  the 
optical  imaging  system  of  the  MNS  be  able  to  locate  the  mine  from  a  safe  distance.  Therefore, 


the  greater  the  stand-off  distance,  the  better  the  vehicle  performance. 

Ln  clear  water  conditions,  where  underwater  visibility  may  cover  quite  large  distances, 
the  system  should  be  able  to  satisfy  the  mission  requirements.  However,  in  many  coastal 
regions  the  underwater  visibility  of  objects  can  be  hampered  by  the  presence  of  absorptive 
and  scattering  matter  in  the  water  column.  In  these  situations,  any  improvement  in  vehicle 
performance  would  be  advantageous,  as  the  chances  for  detonation  of  the  underwater  mine 
would  be  reduced.  In  this  study  we  will  detail  a  suite  of  possible  options  that  ar«  available  to 
the  MNS.  These  options  range  from  some  that  are  extremely  simple  such  as  moving  the  lights, 
to  more  complex  options  which  will  result  in  larger  payoffs  in  improved  vehicle  performance. 

2  Computer  Simulation  of  Underwater  Images 

2.1  The  Physics  of  Underwater  Image  Formation 

In  this  report,  we  will  be  concerned  with  the  propagation  of  electromagnetic  waves 
in  isotropic,  homogeneous  media.  In  this  situation,  the  two  basic  ways  that  light  interacts 
with  matter  are  via  attenuation  and  scattering.  The  attenuation  coefficient  for  pure  water  is 
predicted  via  Maxwell's  equations  to  be  related  to  the  complex  part  of  the  complex  index  of 
refraction.  Absorption  can  also  occur  due  to  the  concentration  of  pigmented  species  in  the 
water.  Scattering  can  be  defined  as  any  divergence  from  a  straight  line  path.  In  the  under¬ 
water  realm  light  can  either  be  diffracted  by  particles  that  are  on  the  order  of  the  wavelength 
of  the  light,  or  refracted  by  particulate  matter  that  has  a  different  index  of  refraction  than 
the  surrounding  sea  or  lake  water. 
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Empirically,  these  phenomena  can  be  grouped  together  via  measurement.  The  total  ab¬ 


sorption  coefficient  c  is  the  decay  constant  that  is  associated  with  the  removal  of  light  intensity 
per  unit  distance: 

Hr)  =  I(o)e~CT  (1) 

Here,  r  is  distance  and  /  is  the  irradiance  at  positions  I(o)  and  I(r). 

The  constant  c  can  be  further  decomposed  as  a  sum  of  two  quantities,  the  absorption  and 
scattering  coefficients  so  that: 

c  =  a  +  b  and  I  =  I(o)e~aT  e~bT  (2) 


The  above  treatment  is  precisely  the  case  for  an  infinitely  small,  thin  beam  where  no  radiation 
is  scattered  back  into  the  beam.  In  this  case,  all  scattering  and  absorption  events  lead  to  the 
loss  of  light  flux.  Transmissometers  measure  an  approximate  value  of  c  by  maintaining  a  very 
small  diameter  beam. 


The  total  scattering  coefficient,  b,  is  the  superposition  of  all  scattering  events  at  all  angles. 
It  can  be  viewed  as  the  integral  of  the  volume  scattering  function  3(6)  over  all  solid  angle: 


b  =  f  3(0)dw  =  27 r  /  3(6)sin0d0 
Jo  Jo 


(3) 


The  four  quantities,  a,  6,  c.  and  3(6),  represent  the  inherent  properties  of  the  medium. 
They  are  distinct  from  the  apparent  properties  as  they  do  not  depend  on  the  radiance  field 
about  a  measurement  point  [Priesendorfer,  1976].  In  principle,  these  quantities  are  all  that 
one  needs  in  order  to  predict  the  propagation  of  light  underwater. 


2.2  The  UNCLES  Computer  Modeling  System 


The  approach  that  we  have  been  pursuing  is  related  to  a  model  formulated  by  McGlamery 
and  coworkers  at  the  Scripps  Institute  of  Oceanography  Visibility  Laboratory  [McGlamery,1979l 
The  model  can  be  described  as  a  hybrid  approach  in  that  linearization  is  incorporated  in  cer¬ 
tain  situations  in  order  to  simplify  the  computation.  When  this  approximation  cannot  be 
made  the  quantities  of  interest  are  calculated  explicitly. 

As  illustrated  in  Figure  2.  an  underwater  imaging  experiment  consists  of  tracing  the 
progression  of  light  from  a  light  source  to  a  camera.  Several  pathways  exist  by  which  the 
light  can  travel  to  the  image  plane  of  the  camera.  These  pathways  may  or  may  not  include 
light  that  has  been  reflected  by  an  object.  Therefore,  the  light  that  has  entered  the  camera 
may  or  may  not  have  been  reflected  by  the  object.  The  light  that  enters  the  camera  without 
reflection  from  the  object  gives  rise  to  an  imaging  component  that  is  called  backscatter. 
Considering  the  light  reflected  by  the  object,  we  can  distinguish  two  components  that  are 
incident  upon  the  camera  plane:  ljlight  that  has  not  been  scattered  in  the  intervening  water, 
called  the  direct  component,  and  2)light  that  has  been  small  angle  scattered,  called  the  forward 
scattered  component.  We  choose  to  represent  an  image  as  being  the  linear  superpostion  of 
three  components:  1  )backscatter,  2)forward  scattering,  and,  3)a  direct  component.  These 
pathways  are  illustrated  in  Figure  2.  Mathematically, 

Ex(total)  =  Ed(direct)  +  E  j ,( f  orwardscatter)  +  Eb,(backscatter)  (4) 
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2.3  The  Direct  and  Forward  Scattered  Component 

In  order  to  calculate  the  light  reflected  from  the  map,  the  irradiance  pattern  incident 
upon  the  reflectance  map  must  first  be  calculated.  The  origin  of  this  irradiance  is  the  light 
source,  considered  here  to  be  a  point  source  characterized  by  a  beam  pattern  ,BP(0.  <^),  a 
function  of  polar  angles  6,  and  p).  In  order  to  calculate  the  radiance  incident  upon  the  map 
this  beam  pattern  is  geometrically  projected.  As  a  final  approximation,  the  scalar  irradiance 
incident  upon  the  reflectance  map  can  then  be  calculated  as: 

e-cR, 

E',(x\y'.  6,  =  BP{6,p>)cosy  — (5) 

Here.  x\y’  refer  to  the  fixed  coordinate  system  with  respect  to  the  planar  reflectance 
map  located  at  z'  =  0.  R„(x,  y,  z,  x\  y\  0)  is  the  distance  from  the  source  to  a  point  on 
the  reflectance  map  and  7  is  the  angle  between  a  perpendicular  to  the  reflectance  map  at 
a  given  x'.y'  location  and  a  line  between  the  x\y'  location  and  the  source.  The  constant  c 
is  the  total  attenuation  coefficient.  The  geometry  is  illustrated  in  Figure  2.  Note  that  the 
primed  coordinate  system  is  associated  with  the  reflectance  map  and  the  unprimed  coordinate 
system  is  associated  with  the  camera  plane.  A  more  accurate  representation  of  the  incident 
irradiancp.  Ej  takes  into  account  the  spreading  of  illumination  due  to  the  small  angle  forward 
scattered  component.  An  approximate  value  can  be  found  via  a  convolution  with  a  point 
spread  function  as  is  typical  in  image  processing  [Rosenfeld  and  Kak:1982].  In  this  case, 

E[(x\  y.  0)  =  E'j(x\  y\  0)  *  g(Ra,G.c,B)  +  E'^x',  y'.O)  (6) 

where 
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g(R,<G,c.  B)  =  [exp{-GR, )  -  exp(-cRs)}f  lexp(-BRsf)  (7) 

and  G  is  an  empirical  constant  (|G|  <  c).  The  operator  F~l  indicates  that  an  inverse 
Fourier  transform  is  taken  of  a  function  B.  an  empirical  damping  factor,  and  f,  a  radial 
frequency  in  cycles/radian.  Wells  [1969]  has  shown  that  in  the  small  angle  scattering  approx¬ 
imation.  the  linear  relationship  of  equation  (3)  is  valid.  This  is  equivalent  to  the  formulation 
proposed  by  McGlamery[1979],  It  is  also  justifiable  on  a  more  theoretical  basis  as  considered 
in  the  treatment  of  Gordon  [1973].  In  Gordon’s  work,  the  relationship  was  found  to  agree 
with  Monte  Carlo  calculations  for  distances  of  up  to  6  attenuation  lengths. 
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In  order  to  compute  the  reflected  radiance,  the  values  of  the  incident  irradiance  pattern 
are  multiplied  by  the  reflectance  values  of  the  reflectance  map.  In  the  UNCLES  system  an 
object  is  represented  by  a  planar  reflectance  map  M(x',y‘ )  located  in  3-dimensional  space 
at  a  given  orientation.  Note  that  M(x',y')  <  1.  Typical  values  for  objects  of  oceanographic 
interest  are  .02  <  M{x',y')  <  .1  [Dixon  et.al.;1983].  In  addition,  a  geometric  factor  of  cost)  is 
taken  into  consideration  in  accordance  with  Lambert’s  law.  The  angle  8  is  the  angle  between 
a  normal  to  the  reflectance  map  at  position  x’,y',and  the  camera  aperture  as  illustrated  in 
Figure  2. 

In  order  to  predict  the  scalar  irradiance  incident  upon  the  image  plane  of  the  camera  it 
is  necessary  to  consider  the  geometric  optics  of  the  camera,  the  attenuation  of  the  medium 
between  the  reflectance  map  and  the  camera,  and  the  spherical  spreading  of  the  reflected 
wave.  Taking  these  factors  into  consideration,  the  scalar  irradiance  incident  upon  the  image 
plane  of  the  camera  can  be  represented  as: 


'i' 

*  S 
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(8) 


r  ■  .  ,  rr  ,  t  t\  ,  r>  \  ^  'U  )  \t i't- r  ^  12 

Ed(x,y)  =  E](x  ,y  )exp(-cRc) - — : - cos  VTt [ — - - j 

4 Jn  R-c 

where  Rc  is  the  distance  from  an  x\y’  position  on  the  reflectance  map  to  the  camera.  /n 
is  the  /  number  of  the  camera  of  focal  length  Fi ,  and  Ti  is  the  transmittance  of  the  lens. 

The  forward  scattered  component  can  then  be  calculated  from  the  direct  component  via 
the  convolution  relationship: 

Efs(x.y)  =  Ed{x,  y.o)  *  g(  Rc.  G,c,  B.)  (9) 

where  g  is  represented  by  equation  [7], 

2.4  The  Backscatter  Component 

The  calculation  of  the  backscatter  component  is  the  most  computationally  intense  part 
of  the  computer  program.  This  is  because  the  linear  approximation  is  no  longer  valid  as 
the  backscattered  light  enters  the  camera  from  a  large  distribution  of  angles.  Alternatively 
stated,  the  small  angle  approximation  used  in  making  the  linear  approximation  in  computing 
the  forward  scattered  component  is  no  longer  valid. 

The  approach  taken  in  this  work  follows  that  of  McGlamery  in  that  3-dimensional  space 
is  sliced  into  planes  of  thickness  Ac'  that  are  parallel  to  the  image  plane  of  the  camera.  The 
irradiance  incident  upon  each  of  these  planes  is  then  calculated,  as  above,  by  computing  the 
direct  component  and  then  adding  an  additional  amount  of  irradiance  due  to  the  light  that 
is  small  angle  scattered. 

The  next  step  in  computing  the  backscatter  component  is  to  determine  the  resulting 
image  due  to  each  of  these  illuminated  slabs.  The  irradiance  incident  upon  the  aperture  of 
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the  camera  is  a  superposition  of  these  illuminated  volume  elements  weighted  by  the  value 


of  the  volume  scattering  function  Representing  the  scalar  irradiance  in  3-dimensional  space 
propagating  away  from  the  light  source  as  E,(x',  y\  r'),  following  the  above  arguments: 


E,(x\  y\  z')  =  E,4{x\  y\ z')  +  Esj,(x\  y\  z') 


Ea,d(x',  y',  z')  =  BP(0.<p) 


exp(—cR3) 


E,j,{x\  y',  z')  =  Ea,d{x\  y\  z’)  *  g(Rs ,  G,  c,  B) 


Now,  taking  a  weighted  superposition  of  the  intensities  from  a  given  slice,  we  can  compute 
the  radiant  intensity  that  is  scattered  toward  the  camera  element  x,y  due  to  volume  AV": 


x ,  y)  -  0(<p)E3(x\  y\  z') AV' 


where  Q(yz)  represents  the  volume  scattering  function  and  AK'  is  an  incremental  volume  in 
3-dimensional  space.  The  angle  <p  is  between  a  line  from  the  volume  under  consideration 
to  the  light  source  and  a  line  from  the  volume  to  the  camera.  Next,  the  image  of  this 
volume  of  water  must  be  calculated.  In  this  case,  following  arguments  similar  to  those  above, 
the  direct  component  of  this  backscattered  scalar  irradiance.  Ebs4{x,y ),  can  be  represented 
[McGlamery,1979]  as: 


Eb3,d(x,y)  =  exP( -cZc,)P(9)E,(x\  y',  z')^^-cos39Ti[Zc'  —  ]2  (14) 

t=l  4*'Tl  ^Ct 

Here.  A Z'  is  the  thickness  of  the  backscattering  volume  AV".  Zcl  is  the  distance  from  a 
point  in  the  camera  to  the  center  of  the  backscatter  slab.  The  value  A’  extends  from  i  =  1 
for  the  first  backscatter  plane  to  N  for  the  plane  adjacent  to  the  target.  As  noted  previously. 


the  total  scalar  irradiance  due  to  the  backscattered  image  ,Eba(x,y),  can  then  be  calculated 


by  adding  an  additional  forward  scattered  component: 

Ebs(x,  y)  =  Eba,d(x i  V)  +  Eba,d{x,y)  *  g(Rc,G,c,B)  (15) 

2.5  Relating  UNCLES  Output  to  System  Performance 

In  evaluating  the  performance  of  an  underwater  imaging  system  there  are  two  important 
criteria  that  must  be  considered.  One  of  these  is  whether  the  amount  of  power  falling  upon 
the  sensing  device  is  large  enough  to  be  able  to  allow  adequate  signal  to  noise  levels  in  the 
image.  In  addition,  the  contrast  in  the  image  must  be  assessed  in  order  to  evaluate  whether 
the  resultant  image  contains  visible  features. 

In  the  case  of  the  power  levels,  adequate  power  incident  upon  the  sensor  can  be  calculated 
via  a  straightforward  conversion  of  UNCLES  output  values.  Several  factors  must  be  taken  into 
account  in  order  to  compute  the  amount  of  photon  flux  that  is  incident  upon  the  aperture  of 
the  camera.  In  this  situation,  we  have  employed  the  approach  taken  by  the  Navy  underwater 
imaging  handbook  [Funk  et.al.  1973]. 

In  order  to  determine  the  power  levels  incident  upon  the  image  plane  of  the  camera,  a 
conversion  formula  has  been  used  to  convert  from  the  output  of  the  UNCLES  program  Mc 
which  is  the  value  for  a  one  watt  lighting  source. 

Ec  =  Mc  *  Ps  *  N,  *  Nw  *  Nr  (16) 

Here,  Ec  is  the  output  value  in  lux,  P,  is  the  amount  of  power  output  by  the  illuminating 
device  integrated  over  1  second  (assuming  a  one  second  exposure),  Ns  is  the  photometric 
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output  of  the  illumination  source,  Nw  is  the  water  radiant  efficiency  transfer,  and  Arr  is  the 
efficiency  of  the  reflector  (only  useful  for  strobe  illumination). 

Standard  values  of  can  be  found  in  the  Navy  handbook.  This  value  is  a  function  of  the 
spectral  characteristics  of  the  lighting  source,  the  distance  between  the  source  and  the  camera, 
and  the  attenuation  constant  of  the  water.  The  value  of  Nw  reflects  how  the  water  passes 
the  different  spectral  bands  of  the  lighting  source.  For  a  blue-green  laser,  the  value  would 
be  larger  than  that  of  a  longer  wavelength  source  which  contains  light  that  is  attenuated 
more  rapidly.  Although  the  UNCLES  system  is  essentially  monochromatic,  by  using  the 
appropriate  values  here,  the  spectral  dependency  of  the  attenuation  can  be  considered  in 
computing  the  incident  power.  Table  1  contains  the  values  of  the  power  conversion  units  that 
should  be  used  for  this  study.  Here  we  have  tabulated  the  conversion  factor  from  UNCLES 
output  to  radiometric  intensity  for  incandescent,  and  thalium  iodide  lighting  sources.  The  left 
hand  column  indicates  the  number  of  attenuation  lengths  and  the  headings  indicate  the  water 
type.  Note  that  the  thalium  iodide  source  has  the  greatest  conversion  efficiency  factor.  This 
is  because  the  spectral  characteristics  of  the  light  source  match  the  pass-band  characteristics 
of  the  water  closely. 

A  more  difficult  area  that  must  be  taken  into  consideration  in  order  to  evaluate  image 
quality  is  the  contrast  transmittance.  Defining  this  quantity  as  the  fraction  of  direct  light 
that  is  incident  upon  the  camera  plane,  it  can  be  represented  as: 


Ed  +  EJs  +  Ebs 


A  minimal  level  of  contrast  is  necessary  in  order  for  the  human  observer  to  discriminate 
target  from  background.  Our  own  experiments  with  the  UNCLES  system  have  led  us  to 


conclude  that  a  value  of  .02  contrast  transmittance  is  adequate  for  viewer  discrimination. 
Furthermore,  an  acceptable  image  will  be  defined  when  this  value  occurs  over  greater  than 
25  percent  of  the  viewing  area.  In  summary,  the  criteria  for  adequate  image  visibility  is:  (1) 
a  power  level  w'hich  allows  the  video  camera  to  image  the  scene  in  its  linear  range,  and  (2)  a 
contrast  value  which  is  greater  than  .02  over  at  least  25  percent  of  the  field  of  view. 

2.6  A  Sample  UNCLES  Run 

In  this  section  a  sample  UNCLES  input  data  set  and  the  resultant  output  of  the  program 
will  be  illustrated.  The  sample  input  data  file  is  contained  in  Table  2.  Here  we  see  that 
the  input  to  the  system  of  computer  programs  consists  of  several  distinct  sections.  The 
first  set  of  input  data  are  the  constants  that  describe  the  environmental  parameters.  The 
total  attenuation  coefficient^  =a  -fb),  the  scattering  to  attenuation  ratio,  the  empirical 
constants  relating  to  the  scattering  power  and  the  loss  of  resolution  as  a  function  of  distance 
are  the  first  four  parameters.  The  volume  scattering  function  (VSF)  is  next  input.  An 
illustration  of  a  slice  through  this  radially  symmetric  function  is  contained  in  Table  2. 

The  next  set  of  input  variables  relate  to  the  characteristics  of  the  camera.  The  x  and  y 
focal  planes,  fn  (f  number  of  the  lens),  and  the  camera  lens  focal  length  are  all  needed  for 
the  computations.  The  beam  pattern  functions  are  input  next.  Here,  a  file  containing  the 
angular  dependence  of  the  power  of  the  beam,  the  total  integrated  power  in  the  beam  (always 
set  to  1  in  the  simulations  considered  here),  the  beam  pattern  width  (used  for  the  striped 
illumination)  and  the  source  pulse  length  in  nanoseconds  for  range  gating  (a  value  of  0  implies 
a  continuous  source)  are  the  next  set  of  variables  that  need  to  be  included. 

In  order  to  specify  the  reflectance  map  pattern  a  2-dimensional  file  is  required.  The  file¬ 
name  is  an  input  variable  to  the  UNCLES  system.  Additional  values  needed  in  order  to 


characterize  the  location  and  the  length  of  the  reflectance  map  are  its  center  x  and  y  coordi¬ 


nates.  its  length  in  the  x  and  y  direction,  and  the  number  of  x  and  y  picture  elements! pixels) 
that  the  map  contains.  Additional  information  is  a  value  for  the  reflectance  of  the  background, 
assumed  to  be  the  area  surrounding  the  reflectance  map. 

Three  sets  of  locations,  and  the  Eulerian  angles  associated  with  the  orientations  of  the 
camera,  and  the  two  sources  are  needed  in  order  to  compute  the  resultant  images.  In  addition, 
the  number  of  rows  and  columns  in  the  final  image  are  input.  Several  options  used  in  the 
computations  also  exist.  Auto-aim  negates  the  need  to  calculate  all  of  the  angles  and  points 
the  lights  and  the  camera  at  the  center  of  the  reflectance  map.  The  options  for  the  convolutions 
are  provided  to  allow  computation  as  Fast  Fourier  Transforms  or  direct  convolutions  in  image 
space.  The  latter  provides  zero  padding  and  eliminates,  any  wrap-around  effects  at  the 
expense  of  computation  time. 

Finally,  a  recent  addition  to  the  program  allows  the  input  of  a  vehicle  which  is  specified 
as  a  set  of  polygons  bounded  by  vertices.  The  vehicle  file  contains  the  number  of  vertices 
and  their  coordinates  in  three  dimensional  space.  This  option  provides  the  ability  to  specify 
baffles  on  the  vehicle  which  can  be  helpful  in  order  to  reduce  the  amount  of  backscatter. 

The  program  is  then  invoked  in  an  interactive,  menu  driven  mode,  and  the  user  is  asked  to 
verify  the  correctness  of  the  input  by  observing  the  predicted  values  of  the  illumination  and 
the  field  of  view  of  the  camera  in  the  target  plane.  If  these  values  are  acceptable  the  program 
is  then  executed  and  an  output  file  containing  matrices,  of  the  illumination  pattern,  the 
direct,  glow,  backscatter  and  total  image  components,  and  the  contrast  transmittance  are  all 
recorded.  We  have  found  that  by  viewing  contour  plots  of  these  matrices  valuable  information 
about  the  distribution  of  the  values  can  be  interpreted.  As  an  example,  considering  the 
simulation  described  above,  a  contour  plot  of  the  total  component  is  illustrated  in  Figure  3 
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and  the  contrast  transmittance  in  Figure  4. 

In  addition,  a  useful  method  of  display  is  to  generate  graphs  of  lines  through  the  two- 
dimensional  data.  For  example,  a  row  or  column  can  be  specified  and  the  values  of  the  total, 
backscatter,  forward  scatter,  and  the  direct  component  can  be  graphed  along  that  coordinate. 
This  is  illustrated  in  Figure  5  for  a  central  row  of  the  image.  The  advantage  of  this  form 
of  presentation  is  that  the  dynamic  range  of  the  image  components  can  be  easily  inferred. 
Finally,  the  program  generates  a  new  parameter  file  which  can  either  be  reused  for  program 
input,  or  saved  as  documentation.  An  additional  note  on  the  parameter  file  is  that  the  scaling 
constants  are  recorded  that  were  used  in  order  to  constrain  the  output  matrix  values  between 
0  and  255. 


3  Vehicle  Characterization  Study 


In  order  to  suggest  changes  to  increase  the  performance  of  the  MNV  system,  the  first 
task  to  be  accomplished  was  to  enter  all  of  the  relevant  system  parameters  into  the  computer 
program.  This  included  the  locations  of  the  existing  lights,  the  cameras, and  the  beam  patterns 
of  the  lights.  This  information  was  obtained  from  the  engineering  drawings  provided  to  WHOI 
by  the  Naval  Coastal  System  Center  (NCSC).  Figure  6  shows  the  configuration  of  the  MNV 
system  as  presented  in  the  MNV  handbook.  Here,  we  can  see  that  two  sets  of  lamps  are 
provided  forward  and  aft.  In  this  study,  we  will  only  be  concerned  with  the  aft  system.  The 
forward  set  of  lamps  are  of  no  interest  for  long  range  viewing  in  that  they  create  a  great  deal 
of  backscatter  which  obscures  the  appearance  of  underwater  objects. 

In  addition  to  the  physical  location  of  the  cameras  and  lights,  the  environmental  circum¬ 
stances  in  which  the  vehicle  might  find  itself  were  considered.  Since  the  vehicle  is  intended 
to  be  used  in  coastal  regions,  the  major  emphasis  was  placed  on  examining  the  system  per- 
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formance  in  poor  visibility  conditions.  Towards  this  end.  a  systematic  study  of  contrast 
transmittance  and  power  requirements  over  the  possible  operating  situations  of  interest  were 
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compiled  using  the  UNCLES  program. 

Figure  7  represents  a  summary  of  the  entire  first  part  of  the  study.  In  this  figure  the 
achievable  altitude  of  the  vehicle  has  been  graphed  versus  the  total  attenuation  coefficient  of 
the  water.  The  points  on  this  graph  were  determined  in  the  following  manner:  first,  a  set  of 
computer  runs  were  performed  at  various  altitudes  for  a  given  water  attenuation  value.  In 
this  case,  we  chose  to  investigate  total  inverse  attenuation  coefficients  of  3.3.  6.6,  and  20  m. 
Then,  the  altitude  at  which  .02  contrast  transmittance  occurred  over  25%  of  the  image  was 
estimated  by  interpolating  the  percentage  values  from  multiple  computer  runs  at  different 
altitudes.  The  results  of  these  computer  runs  are  contained  in  graphs  7a.7b.and  7c.  The  key 
for  all  of  the  simulation  in  this  study  are  listed  in  Table  4. 

In  this  section  we  are  interested  in  the  two  lowest  graphs  in  the  figure.  The  viewing  alti¬ 
tude  is  depicted  on  the  graph  for  the  original  configuration  of  the  vehicle  in  two  situations, 
one  where  the  reflectance  of  the  background  was  .02  and  the  other  where  the  background 
reflectance  value  was  .1.  Information  provided  to  us  from  NCSC  indicated  that  these  val¬ 
ues  would  be  realistic  for  various  underwater  imaging  scenarios.  As  can  be  seen  from  the 
illustrations,  the  practical  performance  limitations  of  the  current  system  for  the  low  contrast 
background  are  to  limit  the  observability  of  the  bottom  to  5  meters  distance  in  3.3  meter 
water,  8.5  meter  distance  in  6.6  meter  water,  and  19  meter  distance  in  20  meter  water.  With 
a  more  reflective  background  of  .1,  the  vehicle  is  predicted  to  be  able  to  obtain  reasonable 
images  at  7  meter  altitude  in  3.3  meter  water,  10  meter  altitude  in  6.6  meter  water,  and  at 
22  meter  altitude  in  20  meter  water.  Thus,  the  practical  limitations  of  the  vehicle  limit  it  to 
1.5-2  attenuation  lengths  in  cloudy  conditions,  and  barely  more  than  one  attenuation  length 


in  clear  water. 


Now  that  the  baseline  configuration  of  the  vehicle  had  been  determined,  a  number  of 
simple  changes  were  considered  whose  aim  were  to  improve  the  performance  of  the  imaging 
system.  The  first  change  that  was  considered  was  to  reposition  the  lights  of  the  vehicle  further 
aft.  A  convenient  place  to  lacate  them  was  found  on  the  vehicle  which  allowed  them  to  be 
placed  approximately  10  feet  back  from  the  camera  in  the  nose  of  the  vehicle.  The  new 
imaging  performance  of  the  vehicle  is  contained  in  Figure  7. 

As  depicted,  a  dramatic  increase  in  the  imaging  performance  was  evident.  In  this  case,  the 
contrast  limited  performance  at  3.3  meters  attenuation  was  now  increased  to  over  7  meters 
for  the  low  resolution  case  and  8.5  meters  for  the  case  of  the  higher  background  reflectance. 
As  can  be  seen  from  the  graph,  in  the  6.6  meter  case  the  new  vehicle  altitude  for  the  low 
contrast  case  is  12.5  meters,  and  in  the  higher  contrast  case,  13.5  meters.  In  the  clear  water 
case,  the  vehicle  altitude  could  be  increased  to  almost  20  meters  in  the  low  contrast  case,  and 
to  29  meters  in  the  high  contrast  case.  Note  that  the  imaging  performance  has  been  increased 
in  all  cases,  however,  substantially  more  in  the  poorer  imaging  conditions. 

The  next  system  modification  that  was  explored  was  to  narrow  the  beam  patterns  for 
the  vehicle  lamps.  The  objective  here  was  twofold:l)concentrate  the  available  lighting  power, 
and  2 (decrease  the  amount  of  backscatter  that  is  incident  on  the  camera  by  decreasing  the 
common  volume  of  intersection  of  the  camera  and  the  lights.  As  can  be  seen  from  the  Figure 
7  system  characterization,  in  all  cases  this  resulted  in  an  increase  in  the  system  performance. 
The  greatest  increase  in  the  imaging  performance  occurred  in  the  cloudiest  water  condition.  In 
the  3.3  meter  water  the  system  performance  increased  to  nearly  9  meters,  a  distance  of  almost 
3  attenuation  lengths  for  the  high  contrast  case.  In  6.6  meter  water  the  imaging  performance 


increased  to  13  meters.  In  the  20  meter  clear  water  case,  the  final  imaging  distance  for  the 
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system  was  37  meters. 


Additional  experiments  were  performed  to  determine  a  set  of  optimal  pointing  angles  for 
the  computer  simulations.  In  this  case,  the  clearest  water  conditions  were  given  priority  in  that 
this  would  result  in  the  longest  range  imaging.  As  can  be  seen  from  the  graph,  this  resulted 
in  maximizing  the  viewing  distance  for  the  20  meter  case,  however,  in  the  cloudier  conditions 
the  change  in  pointing  angles  did  not  result  in  the  highest  image  quality.  If  the  imaging 
system  performance  were  in  fact  optimized  for  the  poorer  water  quality,  some  sacrifice  in 
higher  image  quality  would  result.  Since  the  cloudier  water  conditions  limit  the  performance 
of  the  vehicle  more  severely,  this  is  the  direction  that  would  be  taken  in  a  future  study. 

The  final  study  performed  in  this  first  phase  of  the  research  was  to  examine  the  per¬ 
formance  and  limitations  of  the  lamps  and  video  camera.  Table  2  contains  a  graph  of  the 
required  power  levels  for  the  different  imaging  configurations  in  the  different  water  conditions. 
With  respect  to  this,  we  note  that  the  incandescent  lamps  presently  on  the  vehicle  w-ill  not 
provide  adequate  power  levels  for  some  of  the  extended  range  configurations  that  we  have 
recommended.  Although  adequate  for  the  initial  imaging  range  of  the  vehicle,  according  to 
our  calculations  the  new  extended  configuration  will  be  power  limited.  Two  w'ays  of  combat¬ 
ing  this  power  limitation  would  be  to  increase  the  efficiency  of  the  lamps,  and  to  also  increase 
the  sensitivity  of  the  camera.  We  recommend  that  both  changes  be  implemented. 

In  the  case  of  the  lamps,  thalium  iodide  fighting  sources  provide  two  advantages  for 
underwater  illumination:  One,  they  are  more  efficient  in  that  they  convert  more  electrical 
watts  to  radiant  watts  than  either  incandescent  or  flourescent.  The  second  advantage  is  that 
their  output  spectra  is  much  better  matched  to  the  passband  characteristics  of  the  water. 
These  facts  are  highlighted  in  Table  1,  a  listing  of  conversion  factors  from  UNCLES  units 
which  are  normalized  to  one  watt  output  to  lux.  Here  it  is  seen  that  as  a  function  of  water 
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type  and  distance,  the  thalium  iodide  lamps  can  provide  a  factor  of  5  to  7  increase  in  light 
power  efficency.  In  all  of  the  discussions  that  follow  it  will  be  assumed  that  thalium  iodide 
lamps  of  the  same  wattage  as  the  incandescent  lamps  are  to  be  used. 

The  sensitivity  of  the  existing  camera  has  also  been  examined.  In  Figure  S,  a  graph  of 
signal  output  versus  camera  sensitivity  in  lux,  it  is  noted  that  for  the  newvicon  camera  which 
is  currently  on  the  vehicle,  the  expected  value  of  lux  for  the  extended  configurations  is  almost 
2  orders  of  magnitude  away  from  the  given  sensitivity  of  the  current  earner a( newvicon).  It  is 
therefore  recommended  that  the  camera  be  replaced  with  a  silicon  vidicon  camera  (ultracon). 

In  summary  we  have  shown  that  the  imaging  performance  of  the  vehicle  can  be  approxi¬ 
mately  doubled  in  cloudy  water  conditions  by  observing  the  following  suite  of  suggestions: 

1)  move  the  bottom  lights  aft  to  achieve  maximal  separation  between  camera  and  lights. 
2)  change  the  lights  to  thalium  iodide  with  concurrent  changes  in  beam  patterns(narrowing). 
3-)  change  the  current  camera  to  a  more  sensitive  type( ultracon  is  recommended  for  compat¬ 
ibility). 

4  Advanced  Imaging  Concepts 

4.1  Range-Gated  Imaging 

In  this  section  two  advanced  imaging  concepts  will  be  presented.  As  possible  options  for  the 
MNV  they  have  the  potential  to  allow  longer  range  imaging  than  do  conventional  approaches 
such  as  the  ones  described  above.  Limitations  to  conventional  underwater  lighting  systems 
are  due  to  the  severe  backscatter  component  that  is  present.  The  methods  explored  in  this 
section  have  been  proposed  to  allow  increased  imaging  performance  by  circumventing  the 


large  backscatter  contribution. 


Range-gating  a  light  source  and  a  receiver  is  a  technique  that  has  been  advocated  in  order 
to  d  -rease  the  contribution  due  to  backscatter.  The  basic  idea  is  that  by  propagating  a  short 
pulse  of  light,  a  slab  of  water  of  finite  thickness  will  be  illuminated.  If  the  receiving  device  is 
opened  for  a  time  length  that  is  equal  to  the  length  of  the  light  pulse,  then,  by  varying  the 
onset  time  of  the  device,  information  about  only  a  thickness  of  water  at  a  given  range  will  be 
collected.  This  technique  can  be  used  with  any  pulsed  light  source.  The  idea  is  illustrated  in 
Figure  9. 

We  have  used  an  extended  version  of  the  UNCLES  program  to  examine  the  advantages  and 
disadvantages  associated  with  range-gated  illumination.  The  program  was  extended  in  order 
to  only  allow  backscatter  to  enter  the  camera  from  slabs  of  water  that  were  coincidentally 
illuminated  when  light  was  being  reflected  from  the  target.  Two  sets  of  simulations  were  then 
executed,  one  corresponding  to  a  pulse  length  of  time  equal  to  17  ns  and  another  equal  to  a 
pulse  length  equal  to  35  ns.  In  the  simulations,  it  was  assumed  that  the  receiving  device  could 
be  opened  for  an  exactly  equivalent  amount  of  time.  The  advantage  of  the  longer  pulse  length 
is  that  more  power  can  be  collected  per  pulse  by  the  receiving  device.  The  disadvantage  of 
the  longer  pulse  is  that  the  amount  of  backscatter  that  is  received  by  the  sensor  is  increased. 

An  example  of  the  result  of  using  the  UNCLES  program  to  simulate  the  acquisition  of  a 
range-gated  image  is  shown  in  Figure  10.  In  this  case,  we  have  assumed  coincident  positioning 
of  the  camera  and  the  light  source.  The  inverse  attenuation  coefficient  was  assumed  to  be 
6.6  meters,  and  the  altitude  of  the  device  was  assumed  to  be  28  meters.  The  image  is  for  a 
pulse  length  of  17  nanoseconds.  This  figure  illustrates  the  large  degree  of  image  contrast  that 
can  result  from  the  use  of  a  range-gated  system.  In  fact,  in  all  of  the  cases  considered  in  the 
simulations,  the  contrast  values  were  always  found  to  be  adequate. 


A  set  of  simulations  were  run  for  the  laser  range-gated  case  as  a  function  of  distance  and 


environmental  parameters.  The  attenuation  coefficients  that  were  used  were  identical  to  those 
above.  Power  computations  were  then  performed  in  order  to  determine  the  radiant  intensity 
values  that  would  be  incident  on  the  image  plane  of  the  camera.  For  the  blue-green  laser 
a  power  output  level  of  35mj  was  assumed  with  a  pulse  repetition  rate  of  25  pulses/sec.  In 
this  case,  a  conversion  factor  of  500  watts/lumen  was  used.  Figure  11  contains  a  graph  of 
the  results  of  the  study.  As  can  be  seen,  the  power  levels  for  the  longer  range  images  are  too 
small  to  be  sensed  by  the  receiving  sensor.  Typical  values  for  both  SIT  cameras  and  CCD 
devices  are  5il0_o  lux  [RCA  Handbook;  1974].  This  has  led  us  to  conclude  that  the  current 
possibilities  for  laser  range-gated  imaging  at  large  distances  are  extremely  limited. 

Figure  11  also  contains  the  resultant  power  levels  that  would  be  seen  for  a  gated  strobe 
lamp  with  an  integration  time  of  35ns.  Here,  we  have  assumed  that  the  spectra  of  the  light 
source  would  be  similar  to  that  of  a  Thallium  Iodide  lamp.  It  has  been  assumed  that  strobe 
power  levels  are  a  value  of  250  watts,  and  that  they  are  capable  of  delivering  a  power  of 
85  lumens/watt  of  electrical  power.  In  this  case,  because  the  water  does  not  pass  all  bands 
of  the  broad  spectrum  of  the  light  equally,  Nw  must  be  determined.  This  is  an  empirical 
factor  which  is  a  range  dependent,  environment  dependent  efficiency  factor.  Standard  values 
of  these  coefficients  are  itemized  in  the  Navy  Handbook  . 

Finally,  in  the  case  of  the  strobe  lamp,  the  duty  cycle  for  a  gated,  repetitive,  strobe  needs 
to  be  incorporated  in  the  integrated  illumination  intensity 'calculation.  The  factor  for  a  strobe 
light  repeating  illumination  with  a  pulse  length  of  35  nanoseconds  every  4  microseconds  would 
be  .009.  These  values  are  commensurate  with  using  a  gated  CID  camera[Xybion  Co.;1987]. 
The  main  disadvantage  of  incoherent  gated  illumination,  in  this  respect,  is  that  gating  results 
in  extremely  inefficient  use  of  illumination.  As  can  be  seen,  in  this  case,  the  power  levels  at 
the  receiver  are  less  than  those  due  to  the  laser  light  above.  On  the  basis  of  the  results  of  the 
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computer  simulation  we  note  that  the  laser-gated  and  range-gated  incoherent  illumination 
systems  do  not  hold  the  possibility  for  long  range  viewing.  In  consideration  of  these  results, 
other  imaging  techniques  have  been  formulated  in  order  to  see  whether  improved  performance 
could  result. 

4.2  Scanned  White  Light  Stripe  Performance 

In  this  section  a  new  technique  is  considered  which  consists  of  scanning  a  stripe  of 
incoherent  light  across  the  field  of  view.  The  situation  is  illustrated  in  Figure  12.  Here, 
as  an  example,  the  light  source  and  the  camera  have  been  placed  on  a  vehicle  of  length  5 
meters.  As  is  widely  recognized,  increasing  the  distance  between  source  and  receiver  can 
lead  to  better  images  by  reducing  the  common  backscatter  volume  in  the  final  picture.  The 
proposed  configuration  presented  is  an  extension  of  that  technique.  The  basic  idea  consists  of 
forming  a  narrow  strip  of  illumination  athwartships  and  scanning  this  stripe  along  the  vehicle 
axis. 

The  image  of  the  striped  illumination  pattern  is  then  recorded  on  a  camera  for  a  fixed 
period  of  time.  In  this  example  we  have  used  l/30t/l  of  a  second  as  it  corresponds  to  a  video 
frame  rate.  The  image  can  also  be  acquired  by  an  image  processor  in  the  1  /30</l  of  a  second. 
The  next  step  consists  of  using  the  image  processor  to  detect  the  boundary  of  the  illuminated 
area.  Knowing  the  width  of  the  beam  and  the  approximate  height  of  the  vehicle,  an  area 
of  interest  can  be  delineated  which  corresponds  to  the  area  illuminated  by  the  stripe.  This 
region  of  interest  is  then  numerically  dissected  out  from  the  received  image  and  added  to 
the  contents  of  another  frame  buffer.  Commercially  available  devices  permit  the  feedback  of 
images  into  an  accumulation  image  buffer  at  this  rate. 

This  process  is  then  continued  for  the  number  of  individually  scanned  lines  in  the  picture. 


equal  to  the  total  image  area  divided  by  the  width  of  an  individual  scan  line.  In  this  example. 
20  individual  scan  lines  have  been  delineated  which  are  then  summed.  The  final  image  consists 
of  a  superposition  of  all  of  these  individual  striped  images.  Finally,  contrast  enhancement 


procedures  can  be  used  to  increase  the  visibility  of  such  a  device,  providing  the  receiving 
sensor  has  the  available  dynamic  range.  Figure  12  illustrates  a  flow  chart  for  the  system. 

The  result  of  using  an  extended  version  of  the  UNCLES  program  to  calculate  the  perfor¬ 
mance  of  such  an  imaging  system  is  shown  in  Figure  13.  The  image  simulation  was  performed 
at  an  altitude  of  12  meters  in  3.3m-1  water.  In  this  case,  to  illustrate  the  performance  of 
such  a  device.  20  slices  were  used.  Additionally,  a  gradient  threshold  was  used  to  trigger  the 
delineation  of  the  region  of  interest.  As  can  be  seen,  the  contrast  in  the  image  is  quite  ade¬ 
quate  for  observer  recognition  of  important  features.  In  addition,  the  mean  power  levels  that 
are  obtainable  across  the  object  field  are  adequate  for  camera  detection.  As  illustrated  we 
predict  that  using  light  striped  illumination  image  visibility  to  at  least  4  attenuation  lengths 
is  possible  in  some  situations.  We  feel  that  this  is  a  conser%-ative  estimate.  The  ultimate 
limitations  of  such  a  device  should  be  comparable  to  the  performances  obtained  with  syn¬ 
chronous  scanned  systems,  that  is.  5-7  attenuation  lenrths.  It  is  concluded  that  light  striped 
illumination  has  superior  potential  over  range-gated  techniques  in  certain  long  range  imaging 
situations. 

5  Ambient  Lighting  Study 

At  the  request  of  XCSC  a  study  of  ambient  lighting  situations  that  might  lead  to  useful 
imaging  conditions  was  undertaken.  System  performance  specifications  for  underwater  image 
power  and  visibility  as  a  function  of  water  type  were  requested.  A  modest  amount  of  effort 
was  put  into  a  series  of  simulations  in  the  hope  that  the  UNCLES  mode!  could  be  used 
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to  predict  image  performance  in  various  underwater  imaging  scenarios  where  ambient  light 
would  provide  illumination  for  imaging.  The  results  obtained  were  suggestive,  however  not 
completely  satisfactory  in  that  the  program  was  being  used  in  an  imaging  regime  for  which 
it  was  not  intended.  The  outcome  of  the  simulations  are  presented  here  for  completeness. 

The  geometry  used  to  simulate  ambient  lighting  conditions  was  to  suspend  the  source  at 
an  altitude  of  200  meters  directly  above  the  camera  and  the  target.  In  this  way,  approximately 
parallel  illumination  could  be  simulated.  Although  the  actual  power  in  the  computed  image 
would  not  be  realistic,  it  was  decided  that  the  relative  amount  of  backscatter  could  be  inferred 
and  therefore  that  the  contrast  transmittance  levels  for  the  images  would  be  of  use.  The  results 
of  the  study,  presented  in  Figure  14.  depict  a  range  of  contrast  transmittance  values  contained 
in  the  image  as  a  function  of  water  type  and  altitude. 

The  results  indicate  that  at  an  altitude  of  one  attenuation  length  the  images  will  all 
be  of  high  contrast,  the  contrast  transmittance  values  being  between  .3  and  .5.  When  the 
vehicle  is  at  an  altitude  of  2  attenuation  lengths,  the  3.3  and  6.6  meter  attenuation  constant 
waters  permit  adequate  image  contrast  to  allow  imaging.  In  the  case  of  40  meter  altitude 
and  20  meter  attenuation  coefficient  the  image  contrast  will  be  extremely  low.  The  model 
also  predicts  that  the  images  will  be  of  marginal  contrast  for  the  3.3  and  6.6  meter  imaging 
conditions  at  3  attenuation  lengths  of  altitude. 

These  values  correspond  to  the  case  when  the  sun  is  directly  overhead  and  the  camera  is 
looking  straight  down.  As  a  matter  of  practical  experience,  experienced  divers  have  noted 
that  in  relatively  clear  waters  imaging  distances  can  approach  70  meters.  For  20  meter  water 
conditions  this  corresponds  tc  two  and  a  half  attenuation  lengths.  This  performance  is  greater 
than  that  predicted  by  the  UNCLES  model.  The  most  likely  cause  of  the  discrepancy  is  that 
the  computer  model  is  only  equipped  to  deal  with  singly  scattered  photons.  At  modest 
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water  depths  the  ambient  light  field  is  more  isotropically  distributed  than  that  predicted  by 
the  computer  model.  Less  backscattered  photons  therefore  enter  the  camera  and.  as  a  result, 
more  object  visibility  exists.  A  complete  simulation  of  the  limitations  of  ambient  light  viewing 
would  need  to  take  these  facts  into  consideration.  Under  the  auspices  of  the  present  study 
this  was  not  possible. 

6  Parameter  Sensitivity  Study 

An  important  issue  in  the  examination  of  the  performance  of  underwater  imaging 
systems  is  the  characterization  of  the  environment  in  which  the  light  is  propagating.  As  we 
have  already  seen,  this  requires  knowledge  of  the  total  attenuation  constant,  the  scattering 
to  attenuation  ratios,  the  volume  scattering  function,  and  the  empirical  constants  that  relate 
to  the  spreading  of  the  beam.  A  natural  question  arises  in  this  context  :  What  are  the 
sensitivities  in  the  predicted  images  to  errors  in  the  characterization  of  these  constants? 

In  order  to  examine  the  answer  to  this  question  we  have  performed  a  limited  number 
of  computer  simulations  to  examine  the  variation  in  contrast  transmittance  with  respect  to 
variations  in  the  value  of  the  scattering  to  attenuation  ratio  and  attenuation  constant.  A 
full  multidimensional  simulation  would  necessarily  e  into  account  the  systematic  varia¬ 
tions  in  almost  all  image  parameters  and  the  resultant  changes  in  the  power  and  contrast 
transmittance.  Since  some  of  the  input  parameters  are  complex  this  full  simulation  would  be 
extremely  consumptive  of  computer  power.  Moreover,  a  quick  look  at  some  of  the  physical 
relationships  indicates  which  of  the  parameters  should  be  important. 

The  result  of  a  set  of  simulations  to  characterize  the  change  in  contrast  transmittance 
with  respect  to  scattering  to  attenuation  ratio  are  depicted  in  Figure  15.  Arbitrary  image 
configurations  were  chosen  as  starting  conditions  and  small  perturbations  (10%)  in  the  envi- 


ronmental  parameters  were  performed  in  order  to  examine  the  resultant  changes  in  contrast 
transmittance.  As  can  be  seen  from  the  simulations  the  contrast  transmittance  values  were 
moderately  sensitive  to  changes  in  the  scattering  to  total  attenuation  coefficient  value.  One 
measure  of  this  sensitivity  of  the  change  is  the  slope  of  the  line  graphed  in  Table  3.  Table  3 
outlines  the  slope  versus  water  conditions  for  the  two  sets  of  simulations  performed:  contrast 
transmittance  versus  attenuation,  and  contrast  transmittance  versus  scattering  to  absorbtion 
ratios. 

As  can  be  seen,  the  largest  relative  change  in  contrast  transmittance  to  attenuation  ratio 
was  a  value  of  2.76.  This  infers  that  the  degree  of  accuracy  needed  in  a  for  a  measurement 
of  the  contrast  transmittance  of  accuracy  .01  (half  the  threshold  value)  would  be  .0036.  In 
the  case  of  the  scattering  to  absorbtion  ratio  the  contrast  was  much  less  sensitive,  the  largest 
value  being  approximately  -.7  .  Here,  the  accuracy  needed  in  the  s/a  ratio  would  be  .05 
It  is  therefore  recommended  that  data  collection  devices  be  employed  which  will  yield 
measurements  of  the  total  attenuation  and  absorbtion  to  scattering  ratios  as  described  above. 
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Figure  1 :  Bottom  Mine  Neutralization 
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Figure  7c:  Interpolation  Curves  for  Twenty  Meter  Water 


Figure  8: Camera  Sensitivity  Curves 
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Figure  12:  IMAGE  FORMATION  VIA  STRIPED  ILLUMINATION 
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